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Abstract 

The direct detection rate for supersymmetric cold dark matter (CDM) particles is calculated 
for a number of suitable nuclear targets. Both the coherent and spin contributions are con- 
sidered. By considering representative phenomenologically acceptable input in the restricted 
SUSY parameter space, detectable rates are predicted for some choices of the parameters. The 
modulation effect, due to the Earth's annual motion, has also been considered and found to 
be < 4%. Its precise value depends on the mass of CDM particles (LSP) and the structure of 
the target. 



1 Presented by J.D. Vergados. 



1 Introduction 



There are many arguments supporting the fact that, the cold dark matter of the universe, i.e. 
its component which is composed of particles which were non-relativistic at the time of structure 
formation, is at least 60%. Q There are two interesting cold dark matter candidates: i) Massive 
Compact Halo Objects (MACHO's) and ii) Weak Interacting Massive Particles (WIMP's). The 
MACHO's cannot exceed 40% of the CDM component. |2|, |3[] In the present work we discuss a special 
WIMP candidate connected with the supersymmetry, i.e. the lightest supersymmetric particle 
(LSP). 



We examine the possibility to directly detect the LSP @-[17 via the recoiling of a nucleus (A,Z) 
in the elastic scattering process: 



X + (A,Z) -> X + (Azy 
(x denotes the LSP). In this investigation, we proceed with the following steps: 



(1) 



1) We write down the effective Lagrangian at the elementary particle (quark) level obtained in 
the framework of supersymmetry as described in Refs. f3|j- [|TT 



2) We go from the quark to the nucleon level using an appropriate quark model for the nucleon. 
Special attention in this step is paid to the scalar couplings, which dominate the coherent part 
of the cross section and the isoscalar axial current, which, as we will see, strongly depend on the 
assumed quark model. ]T0| 



3) We compute the relevant nuclear matrix elements |l~8|- ||25|l using as reliable as possible many 
body nuclear wave functions hoping that, by putting as accurate nuclear physics input as possible, 
one will be able to constrain the SUSY parameters as much as possible. 

4) We calculate the modulation of the cross sections due to the earth's revolution around the 
sun by a folding procedure assuming a Maxwell Boltzmann distribution || of velocities for LSP. 



]-0 for LSP detection as e.g. 19 F, 23 Na, 27 Al, 29 Si, m Ca, 
207 Pb has been recently proposed as a theoretical labo- 



There are many popular targets 
73 ' 74 Ge, 127 1, 207 Pb, etc. Among them 
ratory. Furthermore, it can be an important detector, since its spin matrix element, especially the 
isoscalar one, does not exhibit large quenching as that of the light and up to now much studied 29 Si 
and 73 Ge nuclei. 
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Our purpose is to calculate LSP-nucleus scattering cross section using some representative input 

0-0 



in the restricted SUSY parameter space, |TT|-|jT^, |3(], [3T| to compute the coherent LSP-nucleus 
scattering cross sections throughout the periodic table and study the spin matrix elements of 207 Pb, 
since this target, in addition to its experimental qualifications, has the advantage of a rather simple 
nuclear structure. We compare our results to those obtained [[18|] for other proposed cold dark matter 
detection targets. We finally present results obtained by using new input SUSY parameters 
obtained in a phenomenologically allowed parameter space. 
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2 Effective Lagrangian 

Before proceeding with the construction of the effective Lagrangian we will briefly discuss the nature 
of the lightest supersymmetric particle (LSP) focusing on those ingredients which are of interest to 
dark matter. 
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2.1 The nature of LSP 



In currently favorable supergravity models the LSP is a linear combination || |J of the neutral four 
fermions B, W3, Hi and H 2 which are the supersymmetric partners of the gauge bosons B i "" ! 1 1 ' 
and the Higgs scalars Hi and H 2 . Admixtures of s-neutrinos are expected to be negligible. 



n and 



In the above basis the mass-matrix takes the form || |32 

/Mi 
M 2 

-m z cps w m z cpc w 

V m z S{3S w -m z C(3C w -// 



-m z cps w 
m z c p c w 




m z sps w \ 
-m z spc w 

/ 



(2) 



In the above expressions cw = cosOw, sw = sin6\y, Cp = cos(3, s@ = sin/3, where tan/3 = 
(v 2 ) / (vi) is the ratio of the vacuum expectation values of the Higgs scalars H 2 and Hi. u is 
a dimensionful coupling constant which is not specified by the theory (not even its sign). The 
parameters tan/3, Mi, M 2 , \i are determined by the procedure of Refs. 0, |31] using universal masses 
of the GUT scale. 

By diagonalizing the above matrix we obtain a set of eigenvalues rrij and the diagonalizing matrix 
Cu as follows 



/ Br \ 

W 3R 

HlR 

V H 2R j 



{C, 



( XlR \ 

X2R 
X3R 
V X4R J 



( B L \ 

W 2L 
HlL 

V H 2L J 



( XlL \ 

X2L 

X3L 
V X4L ) 



(3) 



Another possibility to express the above results in photino-zino basis 7, Z via 

Ws = sinOwl — cosOy/Z 

B = cosOwl + sin6 w Z (4) 

In the absence of supersymmetry breaking (Mi = M 2 = M and fi = 0) the photino is one of the 
eigenstates with mass M. One of the remaining eigenstates has a zero eigenvalue and is a linear 
combination of Hi and H 2 with mixing angle sin/3. In the presence of SUSY breaking terms the 
B, W3 basis is superior since the lowest eigenstate xi or LSP is primarily B. From our point of view 
the most important parameters are the mass m x of LSP and the mixings Cji,j = 1,2,3,4 which 
yield the xi content of the initial basis states. These parameters which are relevant here are shown 
in Table 1. 

We are now in a position to find the interaction of xi with matter. We distinguish three 
possibilities involving Z-exchange, s-quark exchange and Higgs exchange. 



2.2 The relevant Feynman diagrams 



2.2.1 The Z-exchange contribution 



This can arise from the interaction of Higgsinos with Z which can be read from Eq. C86 of Ref. [3 

■' ^ [HiR^pHin — HiL^/pHiL — (H 2 R r Yi 1 H2R — H 2 l1hH2l)\Z ,j ' (5) 



cosOw 4 
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Using Eq. (^) and the fact that for Majorana particles xipX — 0, we obtain 

L =dU r3l|2 - |C4llW5X1 ^ (6) 

which leads to the effective 4-fermion interaction (see Fig. 1) 

L eff = ^-2(|C 3i r - |C- 41 | 2 )(- 9 - 1 2 Xi7 M 75Xi)-^ (7) 

0056*^4 2cos9 w q z — m z z 

where the extra factor of 2 comes from the Majorana nature of Xi- The neutral hadronic current 
Jf is given by 



J\ = -qi\{\sin 2 9 w - [^(l-7s) - sin 2 e w \r 3 }q (8) 
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at the nucleon level it can be written as 
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Jf = -N lx { sin 2 9 w - g v (^ - sin 2 9 w )r 3 + ^9a1^}N (9) 



Thus we can write 

where 
and 



L eff = --^{xil x l b Xi)Jx{Z) (10) 
J X (Z) = N lx [f v {Z) + f v {Z)r 3 + f° A (ZH + f A {Z) l5 r 3 ]N (11) 



/°(Z) = 2(\C 31 \ 2 -\C 41 \ 2 )^P^- 2 stn 2 6 w (12) 

m z — q z 

f l v {Z) = -2{\C 3l \ 2 -\C A1 \ 2 )-^^g v {\-sin 2 B w ) (13) 

m z — q z 2 



fl(Z) = (14) 

m| 1 
— q z 2 



f\{Z) = 2(|C 31 | 2 -|C 41 | 2 )-^-^ (15) 



with 5v = 1.0, ^ = 1-24. We can easily see that 

fUZ)/fy(Z) = -M^jL- " 1) - -1-15 (16) 

Note that the suppression of this Z-exchange interaction compared to the ordinary neutral current 
interactions arises from the smallness of the mixings C31 and C41, a consequence of the fact that 
the Higgsinos are normally quite a bit heavier than the gauginos. Furthermore, the two Higgsinos 
tend to cancel each other. 



2.2.2 The s-quark mediated interaction 

The other interesting possibility arises from the other two components of xi> namely B and W 3 . 
Their corresponding couplings to s-quarks can be read from the appendix C4 of Ref. They are 

L eff = -gV2{q L [T 3 W 3R -tane w (T 3 -Q)B B ]q L 

- tan8 w q R QB L q R } + HC (17) 



3 



where q are the scalar quarks (SUSY partners of quarks). A summation over all quark flavors is 
understood. Using Eq. (|j) we can write the above equation in the Xi basis. Of interest to us here 
is the part 

L e ff = gV2{{tan9 w (T 3 - Q)C U - T 3 C 2 i)q L XiRQL 

+ tan6 w CnQq R XiLqR\ (18) 

The above interaction is almost diagonal in the quark flavor. There exists, however, mixing between 
the s-quarks q^ and q R (of the same flavor) i.e. 

q L = cosOqqx + sin9qq 2 (19) 

q R = -sinOqqi + cos9qq 2 (20) 



with 



= rnM + ^f) (21) 

m u L ~ m u R , + m tcos2p/ 2 

tan20, = m d (A + ,tan(3) 

d m\ R -m\ + m? z cos2f3/2 v ' 



Thus Eq. ([Tj|) becomes 



L eff = gV2{[B L cos9q q L xm- B R sin9gq R XiL}qi 
+ [B L sin9q-q L XiR + B R cos9q-q R XiL\ fa} 



with 



— —\Ci\tan9u — \C 2 i-, q = u {charge 2/3) 
Bh(q) = — \Cutan9 u) + \C 2 \, q = d (charge —1/3) 
B R (q) = ^tan9 LU Cn, q = u [charge 2/3) 
B R (q) = —\tan9 u) Cu, q = d (charge — 1/3) 

The effective four fermion interaction, Fig. 1, takes the form 

L eff = (gV2) 2 {(B L cos9qq L xiR - B R sin9qq R XiL) 

— z 2 -(B L cos9qX 1R q L - B R sin9qX 1L q R ) 

q l — mqf 

+ (B L sin9 q q L Xi R + cos9qq R Xi L )^—^ — 

q z — mq 2 

(B L sin9qXi R q L + B R cos9qXiLqn)} (23) 

The above effective interaction can be written as 

L eff = L^ff +RR + L L e f f (24) 
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The first term involves quarks of the same chirality and is not much effected by the mixing (provided 
that it is small). The second term involves quarks of opposite chirality and is proportional to the 
s- quark mixing. 



i) The part L L e }+ RR 

Employing a Fierz transformation L^^ RR can be cast in the more convenient form 
L L e fr RR = {gV2?2{-\){\B L \ 2 



, cos 2 9q sin 2 6q _ \ 

— 2 ^2 + — 2 ^)qLl\qLXiRl XiR 

q z — mqf q A — mq 2 



,, sin 2 6 5 cos 2 9 s 

+ \B R 2 ^ + ~2 h )<lRl\<lRXiLl X Xil} 25 

q z — mqf q z — mq 2 



The factor of 2 comes from the majorana nature of LSP and the (-1/2) comes from the Fierz 
transformation. Equation (|25|) can be written more compactly as 

Q 

Left = -^=2{qi\(f3on + P3rt 3 )(1 + 7s)g 



- qix(PoL + P3LT 3 )(l-^)q}(xil X l 5 Xi} (26) 



with 



/4 1 \ 

Por = {qXl R + qXlj \C u tan9 w \ 



\ . i i < i i / r ; ; ; (7 1 1 ' 

PsR = (\xl R -\xl R )\C u tan6 w \ 2 (27) 

1111 

Pol = \-C n tan9 w + -C 2 i\ 2 xl L + \~C n tan0 w - -C 21 \ 2 x\ 

p3L = \-C n tan6w + \c 2 x\ 2 x\ , - \^C n tan8 w - Jc 2 i| 2 x| r 

Z D Z L 



with 



2 2 
2 2 m W , 2 m W 

XqL c <5 9 ~r 5 (5 o 

rriqf-q 2 m^-q 2 

2 2 
2 _ 2 m VK | 2 m 14 / 

Cg = COs9g, Sq = siflQq (28) 

The above parameters are functions of the four-momentum transfer which in our case is negligible. 
Proceeding as in Sec. 2.2.1 we can obtain the effective Lagrangian at the nucleon level as 



L eff +RR = -^(Xil X l 5 Xi)Jx(q) (29) 

Jx(q) = Nix{fv(q) + fv(q)n + flight + fl^h^N (30) 



with 



fv = HPor - Pol), fv = 2(Am - M 
f A = 2g v (p 0R + Pol), f\ = ^9a{Pzr + Psl) (31) 
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We should note that this interaction is more suppressed than the ordinary weak interaction by 
the fact that the masses of the s-quarks are usually larger than that of the gauge boson Z . In the 
limit in which the LSP is a pure bino (Cn = 1, C21 = 0) we obtain 

/4 1 
(3 m = tan 2 e w [-xl R + -xl R 

p 3R = tan 2 9 w (^xl R -\xl R 

tan 2 9 w 2 2 
Pol = 36 (Xu L +Xd L ) 



tan 2 9 w 
36 

Assuming further that Xa R = Xd R = Xa L = Xd L we obtain 



fa = -^(XI L -Xl) (32) 



fv(q)/fv(q) * +1 

tiW/fliq) * +^ (33) 

If, on the other hand, the LSP were the photino (Cn = cos9w, C21 = sinOwiCzi = C41 = 0) 
and the s-quarks were degenerate there would be no coherent contribution (f v = if /3 l = Air)- 

ii) m 

From Eq. (|2"3]) we obtain 



L L * = -( g V2) 2 sin2e^B L (q)B R (q)l-[—^- 
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(qlXirXilQr + cirXilXir<1l) 
Employing a Fierz transformation we can cast it in the form 

As// = + ?7sgXi75Xi - (?<V?) ixi^Xi)) 



where 



+(3-{qr 3 qxiXi + qnihqXHsXi - q° ^qXi^ 1 Xi)] 



(3± = \tan9 w Cii{2sin29u\^Ciitan6 w + \c 2 i]&u 

3 o 2 

T sin29 s [-C n tan9 w - -C 21 ]Aj} 



with 



(m? — m? )M. 2 „ 
U (rn 2 . —q 2 )(m'i —a 2 ) 



■ir 



and an analogous equation for Aj. Here u indicate quarks with charge 2/3 and d quarks with charge 
-1/3. 
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In going to the nucleon level and ignoring the negligible pseudoscalar and tensor components in 



the spirit of Ref. [[33] we obtain 

Leff = %f°(q)NN + fl(q)Nr 3 N]x lX i (34) 

with 

f°(q) = 1.86/3+ and f](q) = 0.48/3_ (Model A) (35) 

(see sect. 2.2.3) The appearance of scalar terms in s-quark exchange has been first noticed in Ref. |7| 
It has also been noticed there that one should consider explicitly the effects of quarks other than u 
and d [ITJ in going from the quark to the nucleon level. We first notice that with the exception of 



t s-quark the q~L — qn. mixing small. Thus 



sin29,r,/\u 



2m u (A + jicotfymyy 
{m\ L - q 2 )(ml R - q 2 ) 

■ no a J 2m d (A + iitan^mly 

sm26 s Ad ~ —- 2 577-5 -77 

(m\ L -g 2 )(m| fl -g 2 ) 

Then the amplitude for this s-quark contribution is proportional to the quark mass (a la Higgs). 
Thus the amplitude for finding such quarks in the nucleon can be computed in a way which is 
similar to that of the Higgs coupling (see Sec. 2.2.3). For the t s-quark the mixing is complete, 
which implies that the amplitude is independent of the top quark mass. Hence in the case of the top 
quark we do not get an extra enhancement in going from the quark to the nucleon level. As we will 
see in the next section we get an enhancement due to quarks other than u and d (see model B in 
the next section). This is not enough, however, to dominate even over /3fy in the SUSY parameter 
space considered here. Thus, fg(q) can be neglected in front of the isoscalar scalar coupling coming 
from Higgs exchange (see sect. 2.2.3). 



2.2.3 The intermediate Higgs contribution 

The coherent scattering can be mediated via the intermediate Higgs particles which survive as 
physical particles (see Fig. 2). The relevant interaction can arise out of the Higgs-Higgsino-gaugino 
interaction which takes the form 

l h X x = (W r H 2 lHI* - W R Hx L Hl* 

- tanB w {B R H 2L Hl* - ~B R H XL Ef)) + H.C. (36) 

Proceeding as above we can express W an B in terms of the appropriate eigenstates and retain the 
LSP to obtain 

L = A^(( C 2i - tan6 w Cn)C 41 XiRXiLH2* 

- (C 21 - tanO^C^XmXiLH?) + H.C. (37) 

We can now proceed further and express the fields H® , ifiT m terms of the physical fields h, H 
and A. The term which contains A will be neglected, since it yields only a pseudoscalar coupling 
which does not lead to coherence. 
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Thus we can write 



where 



with 



dff = -^XXN[fs(H) + ft(H)T 3 ]N (38) 

f°{H) = -(g v + g d )+g a + g c + g b + gt (39) 
fl(H)= l -{g u -g d ) (40) 

9<h = [gi(h)—- + g 2 (H)—-}^, Oi = u,c,t 41 

r ,,.sina /T ^.cosa,m K . , . . 

0* = - 0i fa « + 02 # d— . * = d ' s ' b ( 42 ) 

cosp cosp raw 

gi(h) = A(Ciitan6 w — C 2 i){C4iCosa — Casino) — ^ — ^- (43) 

m z h - q 2 

g2{H) = A{Cntan6 w — C^iC^sina — C^cosa) — ^- — (44) 

m H — q 1 

where is the nucleon mass, and the parameters m/,, m# and a depend on the SUSY parameter 
space (see Table 1). If one ignores quarks other than u and d (model A) and uses m u = 5MeV = 
m d /2 finds [[33 

/° = 1.86(^ + ^/2, fl = 0A9(g u -g d )/2, (45) 
As we have already mentioned, one has to be a bit more careful in handling quarks other than 



u and d since their couplings are proportional to their mass. [|K], [TT]] One encounters in the nucleon 
not only sea quarks (uu, dd and ss) but the heavier quarks also due to QCD effects, which were 
estimated at the one loop level in Ref. [|4], [35] This way one obtains the pseudoscalar Higgs-nucleon 
coupling fs(H) by using effective quark masses as follows 

m u -> f u m N , m d -> f d m N . m s -> f s m N 



m Q 



)W,n, (heavy quarks c,b,t) 



where is the nucleon mass. The isovector contribution is now negligible. The parameters 
/ can be obtained by chiral symmetry breaking terms in relation to phase shift and dispersion 
analysis. [1(], [Uj Following Cheng |34, |35| we obtain 



/„ = 0.021, f d = 0.037, f s = 0.140, Eq fq = 0.240 (model B) 



Another possible solution is 



f u = 0.023, f d = 0.034, f s = 0.400, Eq /q = 0.120 (model C) 



In the present work we will consider these solutions (models B and C) and compare them with the 
solution obtained in the spirit of Addler et al. p3| (Model A above). For a more detailed discussion 



we refer the reader to Refs. fllQ, 11 
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2.3 Expressions for the nuclear matrix elements 



Combining for results of the previous section we can write 

Gf fY.-. - A. 



where 
with 



and 



L eff = -^{(xil X l5Xi)Jx + (XiXi)J} (46) 
Jx = N lx (f° + + /° 7 5 + f\l^)N (47) 



fv = fUZ)+m), fv = fv(Z) + ft® 

f°A = fl(Z)+f A (q), f\ = fl(Z) + fl® (48) 
J = N(f? + fS)N (49) 



We have neglected the uninteresting pseudoscalar and tensor currents. Note that, due to the 
Majorana nature of the LSP, Xi1 X Xi = (identically). We have seen that, the vector and axial 
vector form factors can arise out of Z-exchange and s-quark exchange. @-[f|] They have uncertainties 
in them. Here we consider the three choices in the allowed parameter space of Ref. |3TJ and the 



eight parameter choices of Ref. pi] These involve universal soft breaking masses at the scale. Non- 



universal masses have also recently been employed. | pL0||— [ 1311 111 our choice of the parameters the 



LSP is mostly a gaugino. Thus, the Z- contribution is small. It may become dominant in models 



in which the LSP happens to be primarily a Higgsino. [36] The transition from the quark to the 
nucleon level is pretty straightforward in the case of vector current contribution. We will see later 
that, due to the Majorana nature of the LSP, the contribution of the vector current, which can lead 
to a coherent effect of all nucleons, is suppressed. [Q The vector current is effectively multiplied by 
a factor of /3 = v/c, v is the velocity of LSP (see Tables 2(a), (b)). Thus, the axial current, especially 
in the case of light and medium mass nuclei, cannot be ignored. 

For the isovector axial current one is pretty confident about how to go from the quark to the 
nucleon level. We know from ordinary weak decays that the coupling merely gets renormalized from 
g& = 1 to gjy = 1.24. For the isoscalar axial current the situation is not completely clear. The naive 
quark model (NQM) would give a renormalization parameter of unity (the same as the isovector 
vector current). This point of view has, however, changed in recent years due to the so-called spin 
crisis, [P7|-|]3"9] i.e. the fact that in the EMC data [[JjJ it appears that only a small fraction of the 
proton spin arises from the quarks. Thus, one may have to renormalize f\ by g\ = 0.28, for u and 
d quarks, and g° A = —0.16 for the strange quarks, i.e. a total factor of 0.12. These two 



possibilities, labeled as NQM and EMC, are listed in Tables 2(a), (b). One cannot completely rule 
out the possibility that the actual value maybe anywhere in the above mentioned region. |39| 



The scalar form factors arise out of the Higgs exchange or via s-quark exchange when there is 
mixing [HJ between s-quarks q~i and q~R (the partners of the left-handed and right-handed quarks). 



We have seen in Ref. j|] that they have two types of uncertainties in them. One, which is the most 
important, at the quark level due to the uncertainties in the Higgs sector. The actual values of the 
parameters fg and fg used here, arising mainly from Higgs exchange, were obtained by considering 
1-loop corrections in the Higgs sector. As a result, the lightest Higgs mass is now a bit higher, i.e. 



more massive than the value of the Z-boson. fllq, 17 



The other type of uncertainty is related to the step going from the quark to the nucleon level [ 10 



(see sect. 2.2.3). Such couplings are proportional to the quark masses, and hence sensitive to the 
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small admixtures of qq (q other than u and d) present in the nucleon. Again values of fg and fg in 
the allowed SUSY parameter space are considered (see Tables 2(a),(b)). 

The invariant amplitude in the case of non-relativistic LSP can now be cast in the form 

\M\ 2 = V J 2 + J 2 + J 2 

~ /5 2 |J | 2 + |J| 2 + |^/| 2 (50) 

where m x is the LSP mass, | Jo| and |J| indicate the matrix elements of the time and space compo- 
nents of the current J\ of Eq. (f4"7|) , respectively, and J represents the matrix element of the scalar 
current J of Eq. ([49]). Notice that | Jo| 2 is multiplied by (3 2 (the suppression due to the Majorana 
nature of LSP mentioned above). It is straightforward to show that 

|J | 2 = A 2 |F(q 2 )| 2 (/O-^Z^) 2 (51) 

2 



J 2 = A 2 |F(q 2 )| 2 (fl-fs^^-) (52) 



|J| 2 = 2jiVil< J «ll + ^«i(q)] U)\ 2 (53) 

with -F(q 2 ) the nuclear form factor and 

A A 

n„(q) = E <Kj>~^, «i(q) = E ^(j> 3 (j>^ qXj (54) 

where cr(j), t 3 (j), are the spin, third component of isospin (r 3 |p) = and coordinate of the 
j-th nucleon and q is the momentum transferred to the nucleus. 

The differential cross section in the laboratory frame takes the form [|| 

% = -(— ^TT^^lJoHl-^^^ + lJp+lJl 2 } (55) 

where is the proton mass, rj = m x /m^A, £ = pj • q > (forward scattering) and 

a = —(G F m N ) 2 ~ 0.77 x 10^ 38 cm 2 (56) 
2% 

The momentum transfer q is given by 

q = Qo£, Qo = P— — (57) 
1 + 1] 

Some values of go (forward momentum transfer) for some characteristic values of m x and represen- 
tative nuclear systems (light, medium and heavy) are given in Table 3. It is clear from Eq. fl5"?|) 
that the momentum transfer can be sizable for large m x and heavy nuclei (77 small). 

The total cross section can be cast in the form 

Tfl TV (1 + 1]Y A 

+ (f°M0)) 2 I 00 (u ) + 2/^/in (o)n 1 (o)/oi(«o) 

+ (/]A(0)) 2 /iiN)} (58) 
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The quantities I p entering Eq. (|5"8| ) are denned as (see Fig. 3) 

I p (u ) = (l + p)u (1+p) I"" x 1+ »\F{x)\ 2 dx, p = 0,l 

J 

where F(q) the nuclear form factor and 

u = q 2 b 2 /2 



(59) 



(60) 



Using appropriate expressions for the form factors (in a harmonic oscillator basis with size parameter 
b) we obtain fH-fH 



1 



/pN = T2 { Z 2 I ( z p) z (u ) + 2NZlZ(u ) + N 2 I%(u )} 



r(p) 



r2 r0>) 



where 



^ 'N = (i + p) L ——p 



A=0 u=0 



U 



1+P 



(61) 



(62) 



x 



l-e 



E 

K = 



(a, /3 = N,Z). The coefficients are given in Ref., [|0] for light and medium nuclei, and in 
Ref. j|]J for heavy nuclei. 



The integrals I pp >, with p, p' = 0, 1, (see Fig. 4) result by following the standard procedure of 
the multipole expansion of the e~* q r in Eq. (0). One finds 



I pp '{uq) 



X.K 



O„,(0) ' 



p,p' = 0,1 



(63) 



where, in the special case in which the ground state of 207 Pb is approximated by a 2pi/2 neutron- hole, 
one finds 

loo = hi = hi = 2 I' £ [F 2p (u e)} 2 d£ (64) 
Jo 

Even though the probability of finding a pure 2pi/ 2 neutron hole in the | ground state of 207 Pb 
is greater than 95%, the ground state magnetic moment is quenched due to the 1 + p-h excitation 
involving the spin orbit partners. Hence, we expect a similar suppression of the isovector spin 
matrix elements |43|]-[[|f| (see Table 4(a), (b)). For comparison, we present our results for A = 207 
together with those of A = 29 and A = 73 (Ref. JT§) in Table 4(a). 



3 Convolution of the cross section with the velocity dis- 
tribution 

The cross sections which would be given from an LSP-detector participating in the revolution of 
the earth around the sun would appear retarded. In this section we are going to study this effect 
by using the method of folding. To this aim let us assume that the LSP is moving with velocity v z 
with respect to the detecting apparatus. Then the detection rate for a target with mass m is given 
by 



11 



dN p(0) m 

—r~ — \ \v z \cr{v) (65) 

dt m Am N ] 

where p(0) = 0.3Ge V/cm? is the LSP density in our vicinity. This density has to be consistent with 
the LSP velocity distribution. Such a consistent choice can be a Maxwell distribution || 

fW) = (V^or 3 e-(*>°) 2 (66) 

provided that 

v = ^(2/3)(v 2 ) = 220Km/s (67) 
For our purposes it is convenient to express the above distribution in the laboratory frame, i.e. 

/(v,v £ ) = (^ )-V( v+v -) 2 ^ (68) 

where is the velocity of the earth with respect to the center of the distribution. Choosing a 
coordinate system in which x 2 is the axis of the galaxy, x 3 is along the sun's direction of motion 
(v ) and Xx = x 2 x x 3 , we find that the position of the axis of the ecliptic is determined by the 
angle 7 ~ 29.80 (galactic latitude) and the azimuthal angle uj = 186.3° measured on the galactic 
plane from the X3 axis. |f46| 



Thus, the axis of the ecliptic lies very close to the X2X3 plane and the velocity of the earth is 
ve = v o + v i = v o + v i( sinaxi — cosacosj^ + cosasin'yx.3 ) (69) 

and 

COS (X 

v • vi = v v 1 ft! v vi sin 7 cos a (70) 

V 1 + cor^7 cos z uo 

where Vq is the velocity of the sun around the center of the galaxy, v± is the speed of the earth's 
revolution around the sun, a is the phase of the earth orbital motion, a = 2ir(t — tij/T^, where t\ 
is around second of June and T E = lyear. 



The mean value of the event rate of Eq. (p5|) , is defined by 



Then we can write the counting rate as 

/dN\ _ p(0) m 



\ dt ' m Amjsf 
where 



(^)(E) (72) 



<£>= / ^%/(v,v £ )a(|v|)rf 3 v (73) 



Thus, taking the polar axis in the direction v^, we get 

4 



'0 

or 



67rt>n ^0 J- 1 



v 3 dv / \£\d£e- {v +v z +2vVE V /v °a(v) (74) 



(S) = -jL- r vd vF (^)c-^ + ^a(v) (75) 
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with 



F (x) = X-sinhx - coshx + 1 
One can also write Eq. ([75]) as follows 



2xi /■«> u 



o v 



fx(v)a(v)dv 



with 



7T t>7 



"o 2 



In the case in which the first term in Eq 



becomes dominant, we get 



exp 



[v - v E 



l.'n 



exp 



(v + V E ) 



(76) 

(77) 
(78) 

(79) 



in agreement with Eq. (8.15) of Ref. |3[ In Eq. ( [75] ) the nuclear parameters are implicit in the cross 
section a{y) given from Eq. (|58|). The nuclear physics dependence of (£) could be disentangled 
by taking note of the extra velocity dependence of the coherent vector contribution in a(v) and 
introducing the parameters 



5=^ = 0.27, V--, 
v v 



U = UQlfj 



50) 



where the quantity uq is the one entering the nuclear form factors of Eq. (|59|) for t> = t>o, which in 
this case is given by 

^0 



U 



1 / 2/3 m x c 2 6 \ 2 

2 V (1 + 77) Twy 

Afterwards, we can write Eq. © as 

m x \ 2 (T 



(81) 



TT^N (l+Tj) 2 



A-2Z S 



A 



2?7 + 1 

2(1 + v y 



Jl 



, 2 



+ (/lO (0)) J o + 2^/Afio(0)ni(0)Joi + (/ifii(0)) J n } 



52) 



If we assume that Joo = J01 = «/n, as seems to be the case for 207 Pb, the spin dependent part of 

r n 2 

Eq. ( f78[) is reduced to the familiar expression /^O (0) + f\Qi(0) Jn, where the quantity in the 
bracket represents the spin matrix element at q — 0. 

The parameters J , J p , J pa describe the scalar, vector and spin part of the counting rate, 
respectively, and they are given by 



J (A,W ) 
Jp(A,M ) 



2 e- x ' 2 



6n A 2 Jo 

2 e -A 2 ,00 



67T A 2 Jo 



^e-^Fo(2A^)/o(^ 2 )# 
V 3 e-^F (2A^)/ P (^ 2 )^ 



53) 
54) 
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2 e~^ 2 r°° 2 

/ ^ F (2A^)/ pct (moV 2 )^ 
67T A z Jo 



A = ^ 



1 + Scosasiwy + (5/2) z 



1/2 



35) 
36) 



The parameters I p , I pa have been discussed in the previous section (see Figs. 3 and 4). The 
above integrals are functions of A and u . The latter depends on v , the nuclear parameters and 
the LSP mass. These integrals can only be done numerically. Since, however, A is close to unity, 
we can expand in powers of 5 and make explicit the dependence of these integrals on the earth's 



motion. Thus, 



2 

Jo(A,Mo) = — Bi [ Kq ^ (u ) + 5 sin^y cosa (u ) 
J p (X,u ) = -^=B 2 \Kf\u,)+5sin 1C osaK^ 



37) 



[Uo 



2 

J pa (X,u ) = —=Bx\K^J(u )+5sinjcosaK^J(u ) 

V07T L 



(89) 



The integrals K®, K® and K° pa are normalized so that they become unity at uq = (negligible 
momentum transfer). We find 



Bi 



- I" ^ 2 ''F (2il))dip = -+2v 
e Jo e 



1.860 



Bo 



3e 



3 V e 



Iv) « 4.220 



with 



Furthermore, 



eBi Jo 
2 



K 



3eB 2 Jo 

j poo 



pa 



with F (x) given in Eq. 



eB\ Jo 
|) and 

Fi(x) = 



e~ t2 dt « 0.747 

if)e^Fi(2i/))I (u^ 2 )di(), Z = 0, 1 

V>V^2V>KpW> 2 )#, 1 = 0,1 

#-^(2V0V(^ 2 )#, 1 = 0,1 



— + l)cos/i x - X sin/i x - 1 



The counting rate can thus be cast in the form 



dN 
~dt 



dN 



dt to 



[1 + hcosa) 



(90) 
(91) 

(92) 

(93) 
(94) 
(95) 

(96) 

(97) 



where (^jt)q is the rate obtained from the I = multipole and h the amplitude of the oscillation, 
i.e. the ratio of the component of the multipole I = 1 to that of the multipole I = 0. Below (see 
also Tables 5(a),(b)) we compute separately the amplitude of oscillation for the scalar, vector and 
spin parts of the event rate i.e. the quantity h = Ssiwy K 1 {uq) / K° {uq) . Note the presence of the 
geometric factor siwj = 1/2, which reduces the modulation effect. 
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In order to get some idea of the dependence of the counting rate on the earth's motion, we will 
evaluate the above expressions at w = 0. We get 

K° = K% « K° f = 1 (98) 

R <> = K ~ = TJ^ 53 a402 (99) 

3/(2 e ) + (11 /2W2^ a736 



p 3/e + 7u 



Thus, for siwy ~ 0.5 



Jo « <V = -p= 1.860(1 + 0.054cosa) = 0.857(1 + 0.054cosa) (101) 

V67T 



J p = -= 4.220(1 + 0.099cosa) = 1.944(1 + 0.099cosa) (102) 



We see that, the modulation of the detection rate due to the earth's motion is quite small 
(h « 0.05). The corresponding amplitude of oscillation in the coherent vector contribution, Eq. 
( |101|) , is a bit bigger (h ~ 0.10). However, this contribution is suppressed due to the Majorana 
nature of LSP (through the factor f3 2 ). The modulation due to the Earth's rotation is expected to 
be even smaller. 

The exact K l integrals, for the 1=0 and 1=1, are shown in Figs. 5(a),(b), (c). The most important 
of these integrals, those of Eq. (|93|) associated with the scalar interaction, are shown in Fig. 5(a). 
In Fig. 5(b) we present the integrals of Eq. fl9~4j ) for p = associated with the vector interaction 
(the integral for p = 1 is analogous but it is less important). Finally in Fig. 5(c) the integrals of 
Eq. ( P5| ) for p = 1 and a = 1 are shown. The others are practically indistinguishable from these 
and are not shown (see Ref. 



Before closing this section we should mention that, the folding procedure can also be applied 
in the differential rate in order to obtain the corresponding convoluted expression for da/dQ, i.e. 
before doing the angular integration in Eq. (|55| ) and obtain the total cross section Eq. (0). In 
fact this may be important since the modulation effect in the total cross section is small due to 
cancellations. In fact preliminary results indicate that the modulation effect can get as high as 20%. 
The high value occurs, unfortunately, in the regions where the total differential rate becomes too 
small. 



4 Results and discussion 

The three basic ingredients of our calculation were the input SUSY parameters, a quark model for 
the nucleon and the structure of the nuclei involved. The input SUSY parameters used for the 
results presented in Tables 1 and 2 have been calculated in a phenomenologically allowed parameter 
space (cases #1, #2, #3 of Ref. §0| and cases #4-9 of Ref. [0] 

For the coherent part (scalar and vector) we used realistic nuclear form factors and studied three 
nuclei, representatives of the light, medium and heavy nuclear isotopes (Ca, Ge and Pb). In Tables 
5(a),(b) and 6 we show the results obtained for three different quark models denoted by A (only 
quarks u and d) and B, C (heavy quarks in the nucleon). We see that the results vary substantially 
and are very sensitive to the presence of quarks other than u and d into the nucleon. 
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The spin contribution, arising from the axial current, was computed in the case of 207 Pb system. 
For the isovector axial coupling the transition from the quark to the nucleon level is trivial (a factor 
of qa = 1-25). For the isoscalar axial current we considered two possibilities depending on the 
portion of the nucleon spin which is attributed to the quarks, indicated by EMC and NQM. || The 
ground state wave function of 208 Pb was obtained by diagonalizing the nuclear Hamiltonian 



in a 2h-lp space which is standard for this doubly magic nucleus. The momentum dependence of 
the matrix elements was taken into account and all relevant multipoles were retained (here only 
A = and A = 2). 

In Table 4(a), we compare the spin matrix elements at q = for the most popular targets 
considered for LSP detection 207 Pb, 73 Ge and 29 Si. We see that, even though the spin matrix 
elements Q 2 are even a factor of three smaller than those for 73 Ge obtained in Ref. [18| (see Table 
4(a)), their contribution to the total cross section is almost the same (see Table 4(b)) for LSP masses 
around 100 GeV. Our final results for the quark models (A, B, C, NQM, EMC) are presented in 
Tables 5(a),(b) for SUSY models #l-#3 E0| and Table 6 for SUSY models #4-#9. pT| 



5 Conclusions 

In the present study we found that for heavy LSP and heavy nuclei the results are sensitive to 
the momentum transfer as well as to the LSP mass and other SUSY parameters. From the Tables 
5(a),(b) and 6 we see that, the results are also sensitive to the quark structure of the nucleon. We 
can, however, draw the following general conclusions. 

(i) The coherent scalar (associated with Higgs exchange) for model A (u and d quarks only) is 
comparable to the vector coherent contribution. Both are at present undetectable. For models B 
and C (heavy quarks in the nucleon) the coherent scalar contribution is dominant. Detectable rates 
(dN/dt) > 100 y^Kg^ 1 are possible in a number of models with light LSP. 

(ii) The folding of the total event rate with the velocity distribution provides the total modulation 
effect h. In all cases it is small, less than ±5%. 

(iii) The spin contribution is sensitive to the nuclear structure. It is undetectable if the LSP is 
primarily a gaugino. 
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Table 1. The essential parameters describing the LSP and Higgs. For the definitions see the text. 



Solution 

kJ \J 1 LL U lull 


SI 
ft x 


S2 
ft^ 


S3 
ft° 


SA 

ft^ 


S5 
ft° 


S6 


S7 
ft ' 


S8 


#9 
it 


m (OpV) 
iil x y^jrc v ) 


1 26 


27 


1 02 


80 


124 




1-1 




^0 


m h 


116.0 


110.2 


113.2 


124.0 


121.0 


105.0 


103.0 


92.0 


111.0 


m H 


345.6 


327.0 


326.6 


595.0 


567.0 


501.0 


184.0 


228.0 


234.0 


m A 


345.0 


305.0 


324.0 


594.0 


563.0 


497.0 


179.0 


207.0 


230.0 


tan2a 


0.245 


6.265 


0.525 


0.410 


0.929 


0.935 


0.843 


1.549 


0.612 


tan/? 


10.0 


1.5 


5.0 


5.4 


2.7 


2.7 


5.2 


2.6 


5.3 



Table 2(a). The coupling constants entering C e ff 
#l-#3- 



eqs. (g§), (g7[) and (g9[) of the text, for solutions 



Quantity 


Solution SI 


Solution S2 


Solution S3 


m o 


1.746 x 10 


-5 


2.617 x 10 


-5 


2.864 x 10 


-5 


fvl fv 


-0.153 




-0.113 




-0.251 




P S {H) (model A) 


1.31 x 10- 


5 


1.30 x 10- 


-4 


1.38 x 10- 


5 


fl/f° s (model A) 


-0.275 




-0.107 




-0.246 




fl(H) (model B) 


5.29 x 10" 


-4 


7.84 x 10- 


-3 


6.28 x 10- 


4 


(model C) 


7.57 x 10" 


-4 


7.44 x 10- 


-3 


7.94 x 10- 


4 


fl(NQM) 


0.510 x 10 


-2 


3.55 x 10- 


-2 


0.704 x 10 


-2 


f° A (EMC) 


0.612 x 10 


-3 


0.426 x 10 


-2 


0.844 x 10 


-3 


f\ 


1.55 x 10- 


2 


5.31 x 10- 


2 


3.00 x 10" 


2 



Table 2(b). The same as in Table 1(a) for solutions SA - S9. 



Solution 



#4 



#5 



#6 



#7 



#8 



#9 



(f^fv 
fvl fv 
f°s(A) 
fl(B) 
f!(C) 

f° A (NQM) 
f° A (EMC) 
f\ 



0.225 10~ 4 
-0.0809 
-0.179 10" 4 
-0.531 10~ 2 
-0.315 10~ 2 
-0.207 10~ 5 
6.950 10- 3 
0.834 10- 3 
2.490 10~ 2 



0.190 10~ 4 
-0.0050 
-0.236 10" 4 
-0.145 10~ 2 
-0.134 10~ 2 
-0.407 10~ 5 
5.800 10- 3 
0.696 10~ 3 
1.700 10~ 2 



0.358 10~ 4 
-0.0320 
-0.453 10" 4 
-0.281 lO" 2 
-0.261 lO" 2 
0.116 10^ 4 
1.220 10~ 2 
0.146 10~ 2 
3.440 10~ 2 



0.108 10~ 4 
-0.0538 
-0.266 10" 4 
-0.132 10" 1 
-0.153 10- 1 
0.550 10~ 4 
3.760 10~ 2 
0.451 10~ 2 
2.790 10- 1 



0.694 10" 4 
-0.0464 
-0.210 10-3 
-0.117 10" 1 
-0.118 10" 1 
0.307 10~ 4 
3.410 10~ 2 
0.409 10~ 2 
1.800 10- 1 



0.864 10~ 4 
-0.0369 
-0.131 10-3 
-0.490 10~ 2 
-0.159 10- 2 
0.365 10" 4 
2.360 10- 2 
0.283 10~ 2 
2.100 10" 1 
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Table 3. The quantity go (forward momentum transfer) in units of fm 1 for three values of mi 
and three typical nuclei. In determining go the value (/J 2 ) 1 / 2 = 1CT 3 was employed. 







\ np t \~\ i o 
IN UClcUo 


mi = 30.0 GeV 


mi = 100.0 GeV 


mi = 150.0 GeV 


Ca 


0.174 


0.290 


0.321 


Ge 


0.215 


0.425 


0.494 


Pb 


0.267 


0.685 


0.885 



Table 4(a). Comparison of the static spin matrix elements for three typical nuclei, Pb (present 
calculation) and 73 Ge, 29 Si (see Ref. 0). 



Component 


207 pi, 73 n 29 ca , 
r0 l/2~ tje 9/2+ D< l/2+ 


n?(o) 

fii(0)fi (0) 

ng(o) 


0.231 1.005 0.204 
-0.266 -1.078 -0.202 
0.305 1.157 0.201 



Table 4(b). Ratio of spin contribution ( 207 Pb/ 73 Ge) at the relevant momentum transfer with the 
kinematical factor 1/(1 + rj) 2 , r\ = m/Am N . 



Solution 


#1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


m x (GeV) 


126 


27 


102 


80 


124 


58 


34 


35 


50 


NQM 
EMC 


0.834 
0.645 


0.335 
0.345 


0.589 
0.602 


0.394 
0.499 


0.537 
0.602 


0.365 
0.263 


0.346 
0.341 


0.337 
0.383 


0.417 
0.479 
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Table 5(a). The quantity {dN/dt) in y 1 Kg 1 and the modulation parameter h for the coherent 
vector and scalar contributions in the cases #1 - #3 and for three typical nuclei. 





Vector Contribution 


Scalar Contribution 




(dN/dt) h 


(dN/dt} h 


Case 


(xl0~ 3 ) 


Model A Model B Model C 


#1 
ro ffz 
#3 


0.264 0.029 

n ifi9 n n^Q 
u.iuz u.uoy 

0.895 0.038 


0.151 x 10~ 3 0.220 0.450 -0.002 
n zli n v i n _1 1 49 sfin 1 98 fifin n n9fi 

U.^ilU X 1U l^iZ.oUU IZo.UUU u.uzu 

0.200 x 10- 3 0.377 0.602 -0.001 


#1 
Ge #2 

#3 


0.151 0.043 
0.053 0.057 
0.481 0.045 


0.779 x 10~ 4 0.120 0.245 0.017 
0.146 x 10 _1 51.724 46.580 0.041 
0.101 x 10^ 3 0.198 0.316 0.020 


#1 
Ca #2 

#3 


0.079 0.053 
0.264 0.060 
0.241 0.053 


0.340 x 10~ 4 0.055 0.114 0.037 
0.612 x 10~ 2 22.271 20.056 0.048 
0.435 x 10~ 4 0.090 0.144 0.038 



Table 5(b). The spin contribution in the LSP — 207 Pb scattering for two cases: EMC data and 
NQM Model for solutions #1, #2, #3. 





EMC DATA 


NQM MODEL 


Solution 


(dN/dt) (y-'Kg- 1 ) 


h 


(dN/dt) (y-'Kg- 1 ) h 


#1 


0.285 x 10~ 2 


0.014 


0.137 x 10~ 2 0.015 


#2 


0.041 


0.046 


0.384 x 10~ 2 0.056 


#3 


0.012 


0.016 


0.764 x 10~ 2 0.017 
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Table 6. The same parameters as in Tables 5(a) for Pb for the solutions #4 — #9. Cases #8, #9 
are no-scale models. The values of {dN/dt) for Model A and the Vector part must be multiplied 
by xKT 2 . 







ocaiai 


rait 




Vector 


Part 


Spin Part 






idN\ 
\ dt, /o 




n 


( dNx 
\ n!t /o 


h 




/ dN \ 

\~dT/o 


h 


Case 


A 


D 








EMC 


NQM 




#4 


0.03 


22.9 


8.5 


0.003 


0.04 


0.054 


0.8010 


~ 3 0.16 10- 2 


0.015 


#5 


0.46 


1.8 


1.4 


-0.003 


0.03 


0.053 


0.3710 


^ 3 0.91 10~ 3 


0.014 


#6 


0.16 


5.7 


4.8 


0.007 


0.11 


0.057 


0.4410 


~ 3 0.11 10 2 


0.033 


#7 


4.30 


110.0 


135.0 


0.020 


0.94 


0.065 


0.67 


0.87 


0.055 


#8 


2.90 


73.1 


79.8 


0.020 


0.40 


0.065 


0.22 


0.35 


0.055 


#9 


2.90 


1.6 


1.7 


0.009 


0.95 


0.059 


0.29 


0.37 


0.035 
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FIGURE CAPTIONS 



Fig. 1. Two diagrams which contribute to the elastic scattering of LSP with nuclei: Z-exchange 

in Fig. 1(a) and s-quark exchange in Fig. 1(b). Due to the Majorana nature of LSP only its pseu- 
dovector coupling contributes. J\ can be parametrized in terms of four form factors f v , fy, f%, f\. 
The scalar terms arising from s-quark mixing are negligible in the SUSY parameter space considered 
here. 



Fig. 2. The same as in Fig. 1, except that the intermediate Higgs exchange is considered. This 
leads to an effective scalar interaction with two form factors fg, (isoscalar) and fg (isovector). 



Fig. 3. The integrals I which describe the dominant scalar contribution (coherent part) of the 

total cross section as a function of the LSP mass (m x = mi), for three typical nuclei: Ca, Ge and 
Pb. The value (p 2 ) 1 / 2 = 1(T 3 was used. 



Fig. 4. (a) Plot of the integrals I n as a function of the LSP mass m x = m\. This integral gives 

the spin contribution to the LSP-nucleus total cross section for 207 Pb. The integrals loo an d Ioi are 
similar, (b) Plot of the integrals Iu(u) and Io(u) for Pb. Note that In is quite a bit less retarded 
compared to I . For definitions see the text. 

Fig. 5. Contributions of K integrals (for 1=0 and 1=1) entering the event rate due to earth's 

revolution around the sun: K l in Fig. 5(a), K l in Fig. 5(b) and K l n in Fig. 5(c). The other 
integrals K l 00 and K l 01 are similar to K l n . 
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Fig. 1 



